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Abstract. 

We give the exact expressions of the partial susceptibilities xjj''^ 'S-nd x^f^ 
for the diagonal susceptibility of the Ising model in terms of modular forms 
and Calabi-Yau ODEs, and more specifically, 3F2([l/3, 2/3, 3/2], [1,1]; z) and 
4_F3([l/2, 1/2, 1/2, 1/2], [1, 1, 1];2) hypergeometric functions. By solving the 
connection problems we analytically compute the behavior at all finite singular 
points for xi^^ S'Hd Xd*'- We also give new results for xi^^- We see in 
particular, the emergence of a remarkable order-six operator, which is such that 
its symmetric square has a rational solution. These new exact results indicate 
that the linear differential operators occurring in the n-fold integrals of the Ising 
model are not only "Derived from Geometry" (globally nilpotent), but actually 
correspond to "Special Geometry" (homomorphic to their formal adjoint). This 
raises the question of seeing if these "special geometry" Ising-operators, are 
"special" ones, reducing, in fact systematically, to (selected, fc-balanced, ...) 
qj^\Fq hypergeometric functions, or correspond to the more general solutions of 
Calabi-Yau equations. 



PACS: 05.50.+q, 05.10.-a, 02.30.Hq, 02.30.Gp, 02.40.Xx 

AMS Classification scheme numbers: 341VI55, 47E05, 81Qxx, 32G34, 34Lxx, 
34Mxx, 14Kxx 

Key-words: Lattice Ising model susceptibility, generalized fc-balanced 
hypergeometric fmictions, Saalschiitzian conditions, elliptic curves, elliptic integrals, 
Hauptmoduls, modular equations, modular forms, Calabi-Yau ODEs. 

1. Introduction 



The magnetic susceptibility of the Ising model is defined in terms of the two-point 
spin correlation function as 

oo oo 
M=-oo JV=-oo 

where is the spontaneous magnetization of the Ising model. 

The exact analysis of the Ising model susceptibility is the most challenging and 
important open question in the study of the Ising model today. This study [1] [2] 



began in 1973-76 by means of summing the n"^ particle form factor contribution to 
the correlation function (cco'^a/.tv)- In these papers, it was shown that for T < 

oo 

ksT-xit) = (l-t)'/*- il + J2^(^-)it)), (2) 

n=l 

whereQ t = (smh2E'" /ksT smh2E^ /keT)-^ and for T > T^hy 

OO 

ksT-xit) = {l~t)^^'-Y,x('-+'Ht), (3) 

where t = {smh.2E'" /ksT smh2E^ /ksTf ■ 

The x^"'' are given by n-fold integrals. In [5] the integrals for x^^'' and x*-^' were 
evaluated, and since that time there have been many important studies [31 IHIS] , of the 
behavior as t — >■ 1, of the singularities in the complex i-plane [4l[6l[7] and the analytic 
properties of x^"^ as a function of t for the isotropic case [8]-[T8] for n = 3, 4, 5, 6. 
These studies are still ongoing. 

More recently it was discovered [in] that if in ([Ij the sum is restricted to the 
spins on the diagonal 

oo 

ksT-Xdit) = {KofTw.w) -^'}, (4) 

the diagonal susceptibility reads 

kBT-xa{t)= (^l + f;^?"^w)' (5) 

for T < Tr and 



ksT-xait) = {l~tf"-Y.xT^'\t)^ (6) 



n=0 

for T > Tc. The xi"''Ws 

are n-fold integrals which have a much simpler form than 
the integrals for x*-"' {t) but retain all of the physically interesting properties of these 
integrals. 

For T < Tc, the integrals x'd"\'t) read 



X 



-Q f X2j-l{l- X2j){l''tX2j) 



\ Va;2j(l-a;2i-i)(l -tx2j-i), 

^ n n -tx2j-ix2k)~^ 

l<j<n l<k<n 

^ n {X2]-1 - X2k-l)'^ {X2] - X2k)'^, (7) 
l<j<k<n 

f The classical interaction energy of the Ising model is — j. {E^ i^j^k'^j+i^k + E^'^j,k'^j.k+l) where 
j{k) specifies the row (column) of a square lattice and the sum is over all sites of the lattice. 
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where t is given by t = {sin\Y2E'" /kBTsmh2E^ /kBTy^. 
For T > Tc, the integrals xi^"^^^(0 read 

^n(n+l) /-l /-l 2ji+l 

n+l 

X (^(1 -X2j-i)(l - ta;2j-i) • a;2j-i) 

^ n n -tx2j-ix2ky^ (8) 

l<j<n+l l<fc<n 

^ H (a;2j-i - a;2fe-i)^ {x2j - X2kY ■ 

l<j<k<n+l l<j<k<n 

where t = {smh2E'" /ksT fiin\v2E^ /kBTf ■ For these xi^"^^^ it wiU be convenient, 
in the following, to use the variable X = t^/^ ^ sinh2i;"/fcBTsinh2£''VfcBT. 
In [12] we found that 

xi'^(^) = and: (t) = i . (9) 

and that (^) ^-nd are solutions of differential equations of order 6 and 8. 

The corresponding linear differential operators of each is a direct sum of three factors. 
In both cases, there was a differential equation which was not solved in [19] . 

In this paper, we complete this study of x^^\x) and X^f^ {t) by solving all of 
the differential equations involved. We then use the solutions of these equations to 
analytically compute the singular behavior at all of the finite singular points. In this 
way we are able to give analytic proofs of the results conjectured in Appendix B of [12] 
by numerical means. 

We split the presentation of our results into two parts: the solution of the 
differential equations and the use of the differential equations to compute the behavior 
of X'^\x) and X'^^t) at the singularities. The solution of the differential equations 
is presented in section [5] for x^^"* {x) and in section [3] for x'^'' it) . In particular we 
focus on the difficult problem of solving a particular order-four operator, to discover, 
finally, a surprisingly simple result. The linear differential equation for xi^^(^) is 
studied in section |31 yielding the emergence of a remarkable order-six operator. The 
singular behaviors of x^^'' {x) and xiP (t) are given in section [5] and [6l respectively. 
This analysis requires that the (global) connection problem to be solved. The details 
of these computations are given in appendices C and D. We conclude in section [7] 
with a discussion of the emergence of q+iFq hypergeometric functions, with all these 
previous results underlying modularity [201 121] in the Ising model through elliptic 
integrals, modular forms and Calabi-Yau ODEs 123) . 
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2. Computations for Xd \^) 

It was shown in |19] that xf\x) is annihilated by an order-six hnear differential 
equation. The corresponding linear differential operator £g is a direct sum of 
irreducible linear differential operators (the indices are the orders): 

= Lf)©43)®if . (10) 

The solution of Cl?^ which is analytic at a; = is thus naturally decomposed as a 
sum: 

5o/(£f ) = a['^ ■ xg(x) + ai'^ ■ x%{x) + af • x%{x), (11) 

where the jQ. - are analytic at x = 0. The solutions Xd iC^^) &n.d X^ ^^^) '^^re explicitly 
found in [T^ to be 

xg(x) = and: (12) 

xg(^) = ^^^T^ ■ 2Fi([l/2,-l/2], [1]; x^) - ^ • 2Fi([l/2, 1/2], [1]; x^), (13) 

where one notes the occurrence of x^^^ {x) = xl^lix) in xi^^ (2;)- The last term, x'i^l{x), 
is annihilated by the order-three linear differential operator 



(3) _ n3 , 3 n2{x) 2 _j_ "iW 

no (a;) 



3 - +2d(x) ^ ^ ix + l)ix-l)-x-d{x) " 



where: 



(x-|-l)(a;-l)2.a;.d(x)' 

d{x) = (x -I- 2) (1 + 2 x) (a; + 1) (x - 1) (1 + x + x^) • x, 

no(x) = 2 X® -f 8 x^ - 7 x*^ - 13 x^ - 58 x'' - 88 x^ - 52 x^ - 13 x 5, 

m (x) = 14 x^ -f 71 x^ -f 146 x^ -f 170 x"^ + 38 x^ 

- 112x^ - 94x2 - 19x + 2. 
712 (x) = 8xS + 36x^-|-63x^-H62x3-f 21x2-6x-4. (14) 

(3) 

The linear differential operator has the following regular singular points and 
exponents (z denotes the local variable x — Xs of the expansion around a singular point 
Xg): 

H-x + x2 = 0, p=0, 1,7/2 ^ 

x=0 p = 0, 0, — > ln(z)2 terms, 

x= 1 p=-2, -1, 1^ z"^ (15) 

X = — 1 /3 = 0,0,0 — > ln(z)2 terms, 

X — 00 p=l,l,l — > ln(z)2 terms. 

The last column shows the maximum ln(z)-degree occurring in the formal solutions 

f3) 

of Lg , z being the local variable of the expansion. 
The singularities at x = —2, —1/2 are apparent. 
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By use of the command dsolve in Maple, we found in [T3] that the solution to 
-^3'^'' [Xd-sl ~ which is analytic at a; = is 

- (l^^l)X'+^+l) ■ 3^2([lA2/3,3/2],[l,l];Q), (16) 
where the puUback Q reads: 

Q= ^JSl±lll^ (17) 



Now the coefficients ap'' in the sum decomposition ([TT|) of Xd^''(2;)i can be fixed 
by expanding and matching the rhs of (jlip with the expansion of Xd^"* (2;) ■ This gives 

xHH^) = J-xgl^) -J-xglx). (18) 

By use of a family of identities on 3F2 hypergeometric functions [21] (see eqn. 27 
page 499) the expression (fT6)) of Xd'^3(3^) reduces to 

^S(-> = d'-Ij't-':"'!) ■ [^'''(['/MAl. W; 0)' (19) 

+ ^ ■ 2f'i([l/6,l/31. [1]; 0) - 2Fi([7/6,4/3], |2]; Q) 

It is instructive, however, to discuss further the reason why jQ.-^ix) has this 
solution in terms of 2-F1 functions. 

2.1. Differential algebra .structures and modular forms 

(3) 

From a differential algebra viewpoint, the linear differential operator can be seen 
to be homomorphic[^ to its formal adjoint: 

^3 



^ • adjoint{T2) = T2 ■ adjoint{Lf^), (20) 



where: 



' (l-a;)4 ■ ^ (.T +l)(a; -l)5(2a; + l)(.T + 2)-x ■ " 



and where: 



_ 1 ^^{x) 

4 ' (2a; + l)(a; + 2)(a; +l)(l + x + a;2)(a; - 1)6 -a;' ^ ' 

mi{x) = 2a;6 -6x5 -SSx" -92x3 -81a;2 -34x -6, 

mo{x) = Sx^ - 4x^- 222 a;^ - 769 - 1153a;'' 

- 1341 x^ - 1129 x^ - 490 a: - 84. 

Related to l[20|) is the property that the symmetric squarcQ of actually has a (very 
simple) rational solution R{x). It thus factorises into an (involved) order-five linear 
differential operator and an order-one operator having the rational solution: 

R{x) = Sym^Lf^) = L, ■ [d, - ^ ln(i?(x))). (22) 

§ For the notion of differential operator equivalence see 1251 and 1261 . 

J In general, for an irreducible operator homomorphic to its adjoint, a rational solution occurs for 
the symmetric square (resp. exterior square) of that operator when it is of odd (resp. even) order. 
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In a forthcoming publication we will show that the homomorphism of an operator with 
its adjoint naturally leads to a rational solution for its symmetric square or exterior 
square (according to the order of the operator) . 

Relation (j20p . or the fact that its symmetric square has a rational solution, means 
that this operator is not only a globally nilpotent operator [13], but it corresponds 
to "Special Geometry" . In particular it has a "special" differential Galois group [27] . 
We will come back to this crucial point below, in section [Ql (see (|43|) '). 



(3) 

The operator is in fact homomorphic to the symmetric square of a second 
order linear differential operatoil^ 



X - + i (2x+l)-(x2 + x + 2) 3 1 ^23^ 

' ^ ^ 2 ■ (l + a; + x2).(l + x)-x " + ^ > 



since one has the following simple operator equivalence |25] with two order-one 
intertwinners: 

hf^-Mi = Ni-Sym^{X2), with: (24) 

{l+x)-x , 1 (1 +2x)-{x + 2) 

Ml — —r, TTT- ■ L>x 



[l-xf 2 (l + a; + x2).(l-.T)' 

(l + 2;)-a; 1 24 + 15 a;^ + 8 - 10 a:^ - 69 - 60 a; - 16 

^~ (l-a;)2 ■ 2' (l-a;)3(l +2a;)(x +2)(l + x + a;2) 

The second order operator X2 is not homomorphic to the second order operators 
associated with the complete elliptic integrals of the first or second kind. However, 
from ((20|) and ((22)) . we expect X2 to be "special". This is confirmed by the fact that 
the solution Sol{X2, x) of X2, analytical at a: = has the integrality propertjIJ: 
if one performs a simple rescaling x — > 4 a; the series expansion of this solution has 
integer coefficients: 

Sol{X2,4:x) = 1 +6x2 -24a;^ +60x'* -96a;5 

+ 120 x^ - 672 x'^ + 5238 x^ - 25440 x^ + ■■■ (25) 

From this integrality property [TH |32] , we thus expect the solution of X2 to be 
associated with a modular form, and thus, we expect this solution to be a 2^1 up to 
not just one, but two pullhacks. Finding these pullbacks is a difficult task, except if the 
pullbacks are rational functions. Fortunately we are in this simpler case of rational 
pullbacks, and consequently, we have been able to find the solution [13] to deduce that 
the third order operator L^^^ is 3F2-solvable or 2-^1 -solvable, up to a Hauptmodul [30] 



pullback (see (HI]), ^ 

We can make the modular form character of ()16|) . (|19|) . which is already quite 
clear from the Hauptmodul form of (|17p . very explicit by introducing another rational 
expression, similar to (|17p : 

^ , , 27x • (1 -I- a;) , ^, 

The elimination of x between Q — Q{x) and Qi = Qi{x) (see (fTZl . ([26]) ) gives 
a polynomial relation (with integer coefficients) r((3, Qi) — 0, where the algebraic 
curve r(u, w) = which is, of course, a rational curve, is, in fact a modular curve 

§ Finding X2 (or an operator equivalent to it) can be done by downloading the implementation |28j . 
f See also the concept of "Globally bounded" solutions of linear differential equations by G. 
Christol ESI. 



already encountered |16) with an order-three operator which emerged in (the non- 
diagonal) x'^^ (see [H]): 

-Au^v^ + 12 v"^ ■ {v + u) - 3 u w • (4 -f 4 - 127 u w) 

+ 4:(v + u) ■ (u'^ + v'^ +83uv) - i32uv = 0. (27) 

The hypergeometric functions we encounter in (jl9p . in the expression of the 
solution of ^3 , have actually two possible pullhacks as a consequence of the 
remarkable identity on the same hypergeometric functioij||: 

1 1, 



(1 +2x)- 2i^i([g, 3], [1]; Q{x) 



= {l + x + x^f^- 2F^(\^^-]Al];Qi{x)). (28) 

Other rational parametrizations and puUbacks can also be introduced, as can be seen 
in [Appendix A[ Relation (|28p on 2F1 yields other remarkable relations on the 3i^2 
with the two puUbacks Q (see ([TS]) ) and Qi. their corresponding order-three linear 
differential operators are homomorphic. Consequently one deduces, for instance, that 
3-p2([l/3, 2/3, 3/2], [1,1]; Qi) is equal to the action of the second order operator U2 
on 3F2([l/3,2/3,3/2],[l,l];g): 

{x" +x + lf ■{! - Sx - Sx^)- 3F2([l/3,2/3,3/2],[l,l];Qi) 



3^^2([l/3,2/3,3/2],[l,l];Q) , where: (29) 



U2 



p{x)-[Dl -2-— ln(- — r-^ -A 

^ V ^ dx \{x -l){x +2){l + 2x)J 



p{x) = x^■{\ + xf{\^'2xf{\^%x^\2x^ + %x^^-^x'^), 

q{x) = 8x^° +40x^ +81 a;® -|-84a;^ +24x^ -54a;^ -63x'^ 

- 18a;3 - 2x - 1. (30) 



thus generalizing the simple automorphic relation 

The modularity of these functions can also be seen from the fact that the 
series expansion of (116^ . (llOp . or (1281) have the integrality property |16| . If one 
performs a simple rescaling x Ax, their series expansions actually have integer 
coefficients 1161 1321: 



or 



x'^^liix) = 2 +20 a; +104 -1-560 a;^ -I- 2648 x'^ -j- 12848-^ 



X 



58112 x'' + 267776 a;' + 1181432 a;« + 5281328 a;'' + ••• , (31) 



2Fi([i, [1]; Q(a;))[a; ^ 4-a;] = 1 + 6 a;^ - 24 a;^ + 60 a;^ - 96 a;^ 

+ 120 x^ - 672 a;^ + 5238 a;^ - 25440 a;^ + 81972 a;^° + • • • , (32) 

which can be turned into positive integers if we also change x into — x. 

This provides more examples of the almost quite systematic occurrence [16) in the 
Ising model, of (globally nilpotent [33]) linear differential operators associated with 
elliptic curves, either because one gets straightforwardly elliptic integrals, or because 

I Along this line see for instance | 31| . 
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one gets operators associated with modular forms. For the diagonal susceptibility 
of the Ising model, are we also going to see the emergence of Calabi-Yau-\ike 
operators [HI [53] as already discovered in (the non-diagonal) x*'^'' (see [TB]) ? 

3. Computations for Xd\t) 

We now turn to the computation of Xd'* i^) ? whose differential operator £g'*'' is of order 
eight and is a direct sum of three irreducible differential operators [19; : 

£W = ® L^^) © Lf \ (33) 

The solution of £g'*^ analytic at t ~ 0, is thus naturally decomposed as a sum: 

SoKCt') ^ a['^ ■ x'^lit) + af ■ ^it) + a^' ■ x%it). (34) 

The solutions xjj^j (t) and x|j*2 (0 were explicitly founc^ to be [19] 

xS(^) = and: (35) 

(4) _ 9 (1 +t)-e / 3 5 5 -4t X 

Xd;2W - 8 • (1 _ i)5 • 3^^2([2' 2' 2^' (r~7j^J 



2^^i([l/2, -1/2], [1]; ty - 2Fi([l/2, 1/2], [1]; t) 



(l-t)2 
2+ 

- 2i^i([l/2,l/2],[l];t) • 2Fi([l/2,-l/2],[l];t). (36) 

Here, again, one notes the occurrence of Xd*i (t) which is {t) up to a normalization 

factor. One should be careful that the 3F2 closed form (|36)) for Xd*2(*)' together with 
the previous exact result ([16]), may yield a q+iFq with a rational puUback prejudice 
which has no justification for the moment. 

Similar to the order-three operator for X^^lit), is homomorphic to its 

adjoint and its symmetric square has a simple rational function solution. The exact 
expressions p6l) for Xd*2(^) obtained in a similar way to the solution p^ . ([T^ of 

(3) ' . 

L3 in the previous section. We first find [28^ that the corresponding linear differential 
operator is homomorphic to the symmetric square of a second order operator, which 
turns out to have complete elliptic integral solutions. The emergence in ([36]) of a 
3F2 hypergeometric function with the selectecHI rational pullback —4^/(1 — is 
totally reminiscent (even if it is not exactly of the same form) of Kummers's quadratic 
relation [3TJ [33], and its generalization to 3F2 hypergeometric functions (see the 
relations (4.12), (4.13) in [3S], and (7.1) and (7.4) in [5B^), for example: 



3F2 



([1 + a - /? - 7, 2 ' [1 + « - 1 -}- a - 7]; 



(1 -ty. 

= {1 - tr- 3F2i[a, (3,-f],[l + - 13, 1 + a - ^];t). (37) 

f The first fine in II36I I can, for instance, be found by directly using the command dsolve in Maple, 
and the second line follows by use of identity 520 on page 526 of [24). This result is also easily 
obtained by using Maple to directly compute the homomorphisms between an order-three operator 
and the operator which annihilates 2f'l([l/2, 1/2], 

X The fundamental role played by such specific puUbacks as isogenies of elliptic curves has been 
underlined in |34| . 
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which relates dijferen^ 3F2 hypergeometric functions. In fact, and similar to (j29p . 
we do have an equality between the 3F2 hypergeometric function with the pullback 
u = — 4t/(l — t)'^ and the same 3F2 hypergeometric, where the pullback has been 
changec0 into v — —4(1 —t)/t^, 



3 5 5 -4(l-i) 



3i^2 d, [3, 3]; H ^2 3i^2 b, [3, 3]; , (38) 



3 5 5, „, -4t 



2' 2' 2^' ^ ' ^' (1 -i)2 



where 1^2 is a second-order operator similar to the one in ()29p . The elimination of t 
in these two pullbacks, gives the simple genus zero curve 

u^v^ - 48vu +64- {u+ v) = 0, (39) 

reminiscent of the simplest modular equations |40] . This genus zero curve can also be 
simply parametrized with u — — 4t/(l — t)^ anc ffjl v — 4t • (1 — t). Again, one 
gets an identity, similar to (j38|) . with another order- two intertwinner V2: 



3F2([|,y],[3,3];4t.(l-t))^ (40) 



3 5 5, „, -At 



(l-ty 



3.1. Computation of Xd\{t) 



The third term xi^alO the sum (IMl) is the solution analytic at a: = of the 



!:3 

order-four linear differential operator 



(4) _ ^4 I ^3(0 ^3 I o "2(^) ^2 

4 " * +(t + l)-d4(i) ' (i2-l)-i-d4(t) * 



(i2 - 1) . t . d4(t) * (t- 1) .d4(t)' ^ 

where: 

di{t) = (t^-lOi+l). (i~l)-t, ni(t) = t''- 13t3- 129t2 + 49i -4, 

n2(i) = 5i^-55t^-169t^ + 149t2-28t +2, 
713 (i) = 7t^ -68t^ -llAt^ + 52t -5. 

The operator L^^-* has the following regular singular points and exponents 
i= 0, /9 = 0,0, 0,1 ^ \n{zf terms, 

t— 1, p=— 2,-1,0, 1 — > 2;^^, z^"'^, ln(z) term, 

t = -l, 0,1,2,7 ln(z) term, (42) 

t— 00, p = 0,0, 0,1 In(z)'' terms. 

The last column shows the maximum ln(z)-degree occurring in the formal solutions 
of z being the local variable of the expansion. The singularities at the roots of 

^2 — lOi +1 — are apparent. This order-four operator (|^T|) is actually homomorphic 
to its (formal) adjoint: 

adjoint{L2) ■ L^^'^ — adjoint{L^^^) ■ L2, (43) 

fj Note that the Saalschiitzian difference 15411 (see below) of the 3 F2 at the Ihs of i'STt is independent 
of o, /3, 7 and is equal to 1/2, in contrast with the rhs. 
^ This amounts to changing t into 1 — t. 

ffThis amounts to changing t into — — f) or —1/(1 — t). 
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where L2 is the order-two intertwinner: 

d 



- (A-^ln(Kt))) 



L2 = [Dt - —Mr{t))] ■ Dt where 



.,0 . <''-;°;^t;;i'+" . (44) 

The remarkable equivalence (|331) of operator (HT|) with its adjoint is related to the 
fact that the exterior square of (HT|) has a rational function solution, that is, that this 
exterior square factors into an order-five operator L5 and an order-one operator with 
a rational function solution (which coincides with r{t) in (|44))). 

Exe{Lf) = L,-(Dt- j^Hr{t))). (45) 

In other words, the (irreducible) order-four operator (|4T1) is not only globally nilpotent 
{"Derived from Geometry" [13|), it is a "special" G-operator (33| (Special Geometry): 
its differential Galois group becomes "special" (symplectic or orthogonal groups, see 
for instance p7]). 

This highly selected character of the order-four operator (j4ip is further confirmed 
by the "integrality property" [16] of the series expansion of its analytical solution at 
t = 0: 

5'o/(if') = t + 11/8 -I- 27/16 +2027/1024 1'* 

+ 9269/4096 1^ -t- 83297/32768 1^ + ■■■ (46) 
which, after rescaling t = 16 m, becomes a series with integer coefficients: 

Sol{L^l^) = 16 u +352 li^ +6912 li^ + 129728 li'^ + 2372864 

+ 42648064 + 756609024 + 13286784384 (47) 

+ 231412390144 + 4002962189824 + 68843688570880 

+ 1178125203260416 + 20074611461902336 

+ 340769765322760192^1'' + 5765304623564259328 u^^ 

+ 97249731220784896768^1^ + 1636034439292348588288 + ••• 

This integrality property j32| suggests a modularity |16[ I20[ 121] of this order-four 
operator (|1T|) . The simplest scenario would correspond to pS)) being elliptic integrals 
or, beyond, modular forms that would typically be (up to differential equivalence) 
a 2F1 hypergeometric function with not one, hut two pullbacks (the relation between 
these two pullbacks being a modular curve) . More involved scenarios would correspond 
to Calahi-Yau ODE's [22l [23] and some other mirror maps (see [E]). We have first 
explored the simplest scenarios (elliptic integrals, modular forms), which, as far as 
differential algebra is concerned, amounts to seeing if this order-four operator (|4T|) can 
be reduced, up to differential operator equivalence, to symmetric powers of a second 
order operator. This simple scenario is ruled oui|3- We are now forced to explore the, 
much more complex, Calabi-Yau framework, with two possible scenarios: a general 
Calabi-Yau order- four ODE [53], or a Calabi-Yau order- four ODE that is 4F3 
solvable, the solution like (|46)) being expressed, up to operator equivalence, in term of 
a 4-F3 hypergeometric function up to a pullback that remains to be discovered. This last 
situation would correspond to the 4-^3 Calabi-Yau situation we already encountered in 

t See for instance, van Hoeij's program '28' from ISSAC'2007. 
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X*-^-* (see [H]). The 4F3 solvability is clearly a desirable situation, because everything 
can be much more explicit. 

In constrast with the (globally nilpotent) order-two operators, finding that a given 
order-four operator corresponds to a given 4F3 operator up to a puUback (and up to 
homomorphisms) is an extremely difficult task, because the necessary techniques have 
not yet been developed. Quite often, it goes the other way (no go result): assuming a 
rational pullback, one can rule out a given order-four operator being a ^F^ operator 
with a rational pullback (up to differential operator equivalence). 

In fact, and fortunately, operator (|^T]) turns out to be, a nice example. It has 
singularities at 0, 1, —1, 00, and these points have to be mapped to 0, 1, 00 (i.e. the 
singularities of ^F^ hypergeometric functions), by the pullback. Assuming a rational 
pullback of degree two, there is a systematic algorithm to find all of the rational 
pullbacks of degree two mapping 0, 1, —1, 00 onto 0, 1, 00. This systematic algorithm 
is described in [37] for order-two operators, but the same approach works (with little 
change) for fourth order operators as well. The rational pullback function can actually 
be obtainetfl (with some trial and error) from this mapping of singularities constraint 
and from the exponent-differences, in the same way as in section 2.6 in [39]. The reader 
who is just interested in the surprisingly simple final result and not the mathematical 
structures, in particular the interesting relations between some Calabi-Yau ODEs and 
selected ^F^, can skip the next three subsectiontU and jump directly to the solution 
of gl]) given by ^ with (leT]) . 



3.2. Simplification of L'f^ 

As a "warm up" , let us, for the moment, try to simplify the order-four operator 
(|4ip . getting rid of the apparent singularities — lOi -I- 1 = 0, and trying to take 
into account all the symmetries of (|4ip : for instance, one easily remarks that (j4T|) is 
actually invariant by the involutive symmetry to 1/t. 
Let us introduce the order-four operator 

4 lOx^ -2x-5 „3 , i (5a: + 4).(6x^ -13x +4) ^ 
^ {x-l){l + 2x)x " 4 {x-iy{l + 2x)x^ 

1 .T + 8 3 

^ 4 ■ (a;-l)2(l + 2x)x2 ■ " 4 • (x - 1) (1 + 2 cc) ^3 ' ^ ' 

where 1 + 2 x — is an apparent singularity. One can easily verify that the order-four 
operator (|4ip is the previous operator (|48)). where we have performed the t ^ 1/t 
invariant pullback: 

The operator (|48p is homomorphic to another order-four operator with no 
apparent singularities 

, 5 a; -4 „ 1 (95 - 160 a; + 56) , 
Mi = Dl +2- Dl + - ■ ^ — ■ Dl 

{x-l)- X 4 (a; - 1)2 • a;2 ^ 

1 45a;3 - 124a;2 104 a; - 16 2 a; -5 

^4 (x- 1)3 • x3 " 4- (x- 1)3 . x3' 

% For order-two equations with four singularities (HeunG ...), there are aheady hundreds of cases 
(now all found) , see 1381 . Looking at the size of that table 1381 it is clear that providing an algorithm 
for finding pullbacks will be quite hard. 

t Which correspond, in fact, to the way we originally found the result. 



(49) 
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as can be seen by the (very simple) intertwinning relation: 

M,-D, = (d,+^^^—±^-^).U (51) 
v (a; — 1) (1 + 2 x) • X/ 

This last operator with no apparent singularities, is homomorphic to its adjoint in a 
very simple way: 

adjoint{M4) ■ x'^ ■ {I — x) — x'^ ■ {I — x) ■ M4. (52) 

Do note that, remarkably, the exterior square of Ad^, is an order five operator and 
not the order six operator one could expect generically from an intertwinning relation 
like (|5ip (the exterior square of the order-four operator P5|) is order six with a 
rational function solution (1 + 2x)/x). Taking into account all these last results 
(no apparent singularities, the singularities being the standard 0, 1, c» singularities, 
the intertwinning relation (|5ip . the fact that the exterior square is of order five), the 
order-four operator ([50]) looks like a much simpler operator to study than the original 
operator (|¥T|) . 

3.3. k-balanced 4i^3 hypergeometric function 

Let us make here an important preliminary remark on the ^F^ linear differential 
operators. Let us consider a. hypergeometric function 

4F3([ai, a2, 03, 04], [61, 62, 63]; 0: (53) 

with rational values of the parameters and bj. Its exponents at a; = are 0, 
1 — &i, 1 — 62, 1 — fea, its exponents at a; = 00 are oi, 02, 03, 04, and its exponents 
at a; = 1 are 0, 1, 2 and S where S is the Saalschiitzian difference: 

S ^ (61 + 62 + 63) - (fli + 02 + 03 + 04). (54) 

The Saalschiitzian condition [351 1121 Hll 5 = 1 is thus a condition of confluence of 
two exponents at a; = 1. 

The linear differential order-four operators annihilating the 4-F3 hypergeometric 
function ([55)1 are necessarily globally nilpotent, and they will remain globally nilpotent 
up to pullbacks and up to differential operator equivalenc^ In contrast, the 
corresponding order-four operators are not, for generic (rational) values of the 
parameters and bj, such that they are homomorphic to their (formal) adjoint 
("special geometry"), or such that their exterior square of order-six has a rational 
function solution (a degenerate case corresponding to the exterior square being of 
order five) . 

These last "special geometry" conditions (see (^5)) and (|35|) ). correspond to 
selected algebraic subvarieties in the parameters and bj. In the particular case 
of the exterior square of the order-four operator being of order fiv^ we will show, 
in forthcoming publications, that the parameters ai and bj of the hypergeometric 
functions are necessarily restricted to three sets of algebraic varieties: a codimension- 
three algebraic variety included in the Saalschiitzian condition [331 1131 HI] 5=1 
and two (self-dual for the adjoint) codimension-four algebraic varieties, respectively 
included in the two hyperplanes S = —1 and 5=3. 

X Global nilpotence is preserved by pullbacks (change of variables) and by homomorphisms (operator 
equivalence) . 

t This condition is seen by some authors, see (11) in I41_, as a condition for the ODE to be a Picard- 
Fuchs equation of a Calabi-Yau manifold. These conditions, namely (11) in ^41], are preserved by 
pullbacks, not operator equivalence. 
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Imagine that one is lucky enough to see the order-four operator (which is 
such that its exterior square is of order five) as a 4F3 solvable Calabi-Yau situation: 
one is, thus, exploring particular 4F3 hypergeometric functions corresponding to 
these (narrow sets of) algebraic varieties which single out, particular (A; = —1, 1, 3) 
k-balanced hypergeometric function^ (rather than the well-poised hypergeometric 
functions, or very well-poisec|3 hypergeometric functions j46[ I47j one could have 
imagineqil) . We are actually working up to operator equivalence (which amounts 
to performing derivatives of these hypergeometric functions). It is straightforward to 
see that the n-th derivative of a hypergeometric function shifts the Saalschiitzian 
difference ((54)) by an integer, and that this does not preserve the condition for 
the exterior square of the corresponding order-four operator to be order five: 
it becomes an order-six operator (homomorphic to its formal adjoint) with a 
rational function solution. The natural framework for seeking ^F^ hypergeometric 
functions (if any) for our order- four operators dH]), pSj) . ([5D|) is thus (selected) k- 
balanced hypergeometric functions (rather than the well-poised, or very well-poised, 
hypergeometric functions |46] ...). 

3.4- L^i^ is iFy, solvable, up to a pullhack 

Let us restrict ourselves to the, at first sight simpler, order-four operators ([50)) : 
even if we know exactly the rational values of the parameters a; and bj , finding that 
a given order-four operator corresponds to this given ^Fy, operator, up to a puUback 
(and up to homomorphisms), remains a quite difficult task. We have first studied the 
case where the puUback in our selected 4-^3 hypergeometric functions is a rational 
function. This first scenario has been ruled out on arguments based on the matching 
of the singularities and of the exponents of the singularities. 

We thus need to start exploring puUbacks that are algebraic functions. Algebraic 
functions can branch at certain points (this can, for instance, turn a regular point into 
a singular point). The set of algebraic functions is a very large one, so we startecU 
with the simplest algebraic functions situation namely, square roots singularities. A 
first examination of the matching of the singularities, and of the exponents of the 
singularities, indicates that we should have square roots at a; = 1 only. 

Along this square root line, let us recall the well-known inverse Landen 
transformation in terms of fc, the modulus of the elliptic functions parametrizing the 
Ising model: 



% fc-balanced hypergeometric functions correspond to the Saalschiitzian difference I I54I I being an 
Integer k: S = k. 

tt Note that very well-poised hypergeometric series are known |45l to be related with f(2),f(3), 
which are constants known to occur in the Ising model 

f Note that the conditions to be well-poised hypergeometric series are actually preserved by the 
transformation — > 1 — a; , fej — > 2 — bj, which corresponds to changing the linear differential 
operator, associated with hypergeometric functions, into its formal adjoint. 
X And also because we had a Ising model prejudice in favour of square roots (18| ... 



k 



1 - ^/T^ 



(55) 



1 -t- VI - A:2' 
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In terms of the variable x — fc^, this inverse Landen transformation reads: 



VI + vr^/ 



- -4-(2-x)- ^ J . (56) 



^ -7^^ = ^ = (58) 



Using this puUback P{x), we have actuahy been able to obtain the solution of the 
order-four differential operator (|50p in terms of four terms like 

,F,{[^,^,^,^],[l,l,l];P{x)). (57) 

This slightly involved solution is given in [Appendix B[ 

We can now get the solution of (pij) . the original operator L'f\ from this slightly 
involved result, since (|4T|) is (|50|) up to a simple puUback namely the change of variable 
(|49p. Coming back with (j49|) . to the original variable i in L^f\ the previous puUback 
(|56|) simplifies remarkably: 

At \ _ 1 + 1 - 

(1 -i)2j ^ 2-t2 2-t2 

the Galois conjugate of (|56l) giving l/t^. Of course, once this key result is known, 
namely that a puUback works, it is easy to justify, a posteriori, this simple monomial 
result: after all, l'"^^ has singularities at 0, 1, —1, oo, and these points can be mapped 
(under t^) to 0, 1, oo (i.e. the singularities of 4F3 hypergeometric functions). 

Fullbacks have a natural structure with respect to composition of function^ It is 
worth noting that (j58p describes the composition of two well-known isogenics of elliptic 
curves, the inverse Landen transformation (|56p . and the rational isogeny —4t/(l — t)'^ 
underlined by R. Vidunas [3X^ and in ^34) , giving the simple quadratic transformation 
t t^. 

All this means that the solution of L^^^ can be expressed in terms of the 
hypergeometric function 

4^3([^, ^, ^, [1, 1, (59) 

and its derivatives. Actually considering the hypergeometric operator % having ((5^1 
as a solution, it can be seen to be homomorphic to (|41l) 

A^-n - if -Ag, (60) 

where the order-three intertwinners A3 and read, respectively, (with d'i{t) = 
t-{t + l)-{t~lf - {t^ -lQt + 1)): 

A3 = 2-{l+t)-t^ -Dl +- ^-^ t^-Dl 

1 31*2 _ 4t _ 11 
+ 3 t-Dt +t, (61) 

t Suppose that an operator O2 is a puUback of an operator Oi, where the puUback / is a rational 
function and that O3 is also a puUback of Oi , where the puUback is a rational function g. Then O3 
is also a puUback of 02- To compute this puUback function, one has to compose g and the inverse 
of/. 
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3 t-d-iit) 
1 1 



A3 = --Df +-■——■ {lOt^ -107 ~ 225t^ + 163t~ 17) ■ 

t — 1 3 ci^yt) 

^ ^ {5t^ ~ 66t^ - 900f + 290t- 33)- Dt 

, 3 (.3-21t^+99t- 23). (62) 

From the intertwinning relation (|60| . one easily finds that the sohition of L^^"^ 
which is analytic at t = 0, is A3 acting on (|59|) : 

= A3[4F3([l/2, 1/2, 1/2, 1/2],[1, 1, l];i2)j. (63) 

Having with (j63l) a normalization for xl^\, we can now fix the values of the 
coefficients in the sum (|34l) for x|;*''(^)- They can be fixed by expanding and 
matching the rhs of (l34l) and xli\t)- This gives 

xfHt) = J3.xS(t) +^-xSW -^-xSW- (64) 

Remark: It is quite surprising to find exactly the same 4F3 hypergeometric 
function (|59p with the exact same remarkably simple puUback t^, as the one we 
already found [16 in the order-four Calabi-Yau operator L4 in x^^\ 

Comment: Of course, from a mathematical viewpoint, when looking for a 
puUback, one can in principle always ignore all apparent singularities. These 
calculations displayed here look a bit paradoxical: the calculations performed with 
the (no apparent singularities) operator (jSOp . which looks simpler (it has an exterior 
square of order five, and is very simply homomorphic to its adjoint, ...) turns out to 
have a more complicated puUback (j56p . than the amazingly simple puUback (namely 
t^) wc finally discover for the original operator (|4T]) (see (1631) ). The "complexity" of 
the original operator (|4ip is mostly encapsulated in the order-three intertwinner ^43 
(see ([6T|) ). The "4F3-solving" of the operator amounts to reducing the operator, up to 
operator equivalence ([60]) . to a 4F3 hypergeometric operator up to a puUback. Finding 
the puUback is the difficult step: as far as "4i^3-solving" of an operator is concerned, 
what matters is the complexity of the puUback, not the complexity of the operator 
equivalence. 

Ansatz: Of course, knowing the key ingredient in the final result (1551) . namely 
that the puUback is just t^, it would have been much easier to get this result. Along 
this line one may recall the conjectured existence of a natural boundary at unit circle 
\t\ — 1 for the full susceptibility of the Ising model, and, more specifically for the 
diagonal susceptibility n-fold integrals we study here, the fact that the singularities 
are all iV-th roots of unity (TV integer). Consequently, one may have, for the Ising 
model, a t'^ {N integer) prejudice for puUbacks. 

In forthcoming studies of linear differential operators occurring in the next (bulk) 
x(")'s or (diagonal) X^^^'s, when trying to see if these new (Calabi-Yau like, special 
geometry) operators are q+iFq reducible up to a puUback, we may save some large 
amount of work by assuming that the corresponding puUbacks are of the simple form 
where N is an integer. 
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("5) 

4. The linear differential equation of ' in mod. prime and exact 
arithmetics 

Tlie first terms of the series expansion of Xd'^ i^) read: 

v(5)f^^ - ^ ^.12 I 39 14 5085 ig 

^'^ ^ ' 262144 ^ 1048576 ^ 67108864 ^ ^ 

9 17 33405 lo 315 i„ 

-I . X -\ • X -\ • X + • • • 

67108864 268435456 536870912 

where x = t^^"^ — svah2E^/kT s'mh2Eh/ kT is our independent variable. In order 
to obtain the linear differential equation for Xd\^)j have used in |19j a "mod. 
prime" calculation which amounts to generating large series modulo a given prime, 
and then deduce, the linear differential operator for xj^^-* (x) modulo that prime. With 

3000 coefficients for the series expansion of xli^ i^) modulo a prime, we have obtained 
linear differential equations of order 25, 26, • • • . The smallest order we have reached is 
19, and we have assumed that the linear differential equation of Xd\^) is of minimal 
order 19. In [18j . we have introduced a method to obtain the minimal order of the 
ODE by producing some (> 4) non minimal order ODEs and then using the "ODE 
formula" (see [Tl[15l [18] for details and how to read the ODE formula). The ODE 
formula for xf"^ i^) reads 

31Q + 191? - 302 = (g + l)-(i:> + l) (66) 

confirming that the minimal order of the ODE for x^^^ (x) is 19. Note that the degree of 
the polynomial carrying apparent singularities should be 237 (see Appendix B in [Ml). 
Call C^iQ the differential operator (known mod. prime) for xli^^)- The singularities 
and local exponents of C^g ar^ 

x^O, p= 05,1/2,1^,23,43,3,7,12, 

a; = oo, P = l^ 3/2, 2^ 3^, 4, 5^, 8, 13, 

x^l, p= -3,-2,-l,0^2^4^ ••• , (67) 

a; = -i, p= o^2'^,4^6^8^lo^ ••• , 

x^xo, p= 5/2,7/2,7/2, ••• , 

X = xi, p — 23/2, • • • 

where xq (resp. xi) is any root of 1 + 2; + = (resp. 1 + x + x"^ + x^ + x^ = 0), 

and the trailing • • • denotes integers not in the list. Note that, in practice, we do 

(5) 

not deal with the minimal order differential operator but with an operator of 
order 30 (that c'fg rightdivides) : order 30 is what we have called in [Ml [151 EH] 
the "optimal order" , namely the order for which finding the differential operator 
annihilating the series requires the minimum number of terms in the series. With 
the tools and methods developed in [M] [151 [l2]i we are now able to factorize the 
differential operator and recognize some factors in exact arithmetic. This way, we 

(5) 

may see whether some factors occurring C\q follow the "special geometry" line we 
encountered for x'"f''ix) and X^f\t). Our first step in the factorization of C^^g is to 
check whether c'^^ (the differential operator for xf^) is a right factor of c'fg , meaning 
that the solutions of £g'^^ (and in particular the integral xliH^)) ^^re also solution of 

J The local exponents are given as (e.g.) 2^ meaning 2, 2, 2. 
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rfg) = Lf^®L^^ (69) 



C^iQ ■ This is indeed the case. Using the methods developed in [HI [HI [IT], we find 
that the series for the difference xf^ (x) — a x'f' {x) requires an ODE of minimal order 
17 for the valu^ a — S. This confirms that is in direct smri in cf^, and that 
some (order-four) factors of are still in dy): 

£fg) - Cf^ ® £^'7^ (68) 

This order four factor is obviously L\' ® L\' . Since these factors are in direct sum 

in £g'^'' , the order-seventeen operator is also the annihilator of x'^ i^) ~ P X'd^ (^) 
for /? = 4, meaning that we also have 

-19 

At this step, the differential operator is known in prime. To go further in the 
factorization, we use the method developed in [141 115) along various singularities and 
local exponents of cf.^ which rea4i|: 

a; = 0, /9= 0^1/2,1^2^4^3,7, Mzf,z^l'^, 

x^oo, p= 1^3/2,2^3^4,5^8, In(z)'^, 2^/2^ 

x^l, /9= -3,-2,-1, 0"^, 23,- •• , ln(z)^ 

a; = -l, p= 0^2■*,4^6^8^••• , ln(z)^ (70) 

x = xo, p= 5/2,7/2,7/2,- •• , z^'"^ , 1'^ , z'^ '"^ \n{z), 

x^xi, p= 23/2,- •• , z^^/^ 

where xq and xi are again the roots of l + x + x^ = and 1 + x + + x'^ + x'' = 0, and 
the trailing • • • denotes integers not in the list. The last column shows the maximum 
ln(z)-degree occurring in the formal solutions of , z being the local variable of the 
expansion. 

Use is made of section 5 of |19j to recognize exactly some factors. This is 
completed by an usual rational reconstruction |17j . We are now able to give new 
results completing what was given in section 5 of T9| . The linear differential operator 
dy) has the factorization: 

- L'i'-C[\\ (71) 

The linear differential operator Cl^ has been fully factorized and the factors are 
known in exact arithmetic (the indices are the orders) 

£(5) = Lf) ® e (w^f ) . [/f )) e (if • v}'^ ■ t/f )), (72) 

and are given in [Appendix C| The factor if is the only one which is known in 
prime^ and it is irreducible. The irreducibility has been proven with the method 
presented in section 4 of [M]. This is technically tractable since there are only two 
free coefficients (see ([77)) below) that survive in the expansion of the analytical series 
at X = of if. In the factorization ([7T|) . ([7^ of £f and Cy^, the factors are 

f3l f3l 

either known and occur elsewhere {L\ , ) 01 simple order-one linear differential 

^ Comparing with eq. (58) in |19| . one should not expect a 1/2 contribution, since the sum on the 
g^{N,t)'s still contains xf . 

f There are solutions analytic at a; = with exponents 0, 1, 2, 4, 7. 

(5) 

t This factor known only in primes does not allow the computation of the singular behavior of xjj ■ 
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operators {u'l^\ V^/^'', M^i^'*), except the order-four operator L^f^ and the order-six 
operator ig^''. It is the 
are "Special Geometry" 



(5) 

operator ig . It is then for these specific operators that we examine whether they 



4-1. The linear differential operator L^^'^ 

The order-four finear differential operator L\ has the foUowing focal exponents 

a; = 0, p --2,-2,-1,0, 

X = oo, p = 3, 3, 4, 5, 

x^l, p = -2,-2,-2,-2, (73) 

X = -I, p = -2,-2, 0, 0, 

1 +x +x^ = 0, p==-l,0,l,2. 

At all these singularities xq, the solutions have the maximum allowed degree of log's 
(i.e. ln(a; — xqY), except at the singularities roots of 1 + x + a;^ = 0, where the 
solutions carry no log's. In view of the negative local exponents in (j73p . we introduce: 

Pl[x) ^ (S • (a; - 1)2 . (a; + 1)2 . (.t2 + x + 1)) \ (74) 

Then, if we consider the linear differential operator /i(a;)^^ • L^^'^ ■ fi{x), nothing 
prevents (as far as the p's and log's are concerned) to check whether this conjugated 
operator is homomorphic to a symmetric cube of the order-two linear differential 

(5) 

operator of an elliptic integral. We find the solution of L\' , which is analytic at 
X — 0, as a cubic expression of 2F1 hypergeometric functions, with palindromic 
polynomials 



Sol{Lf^) = Ai(a;) • a;^ • + + 8 • (2 + 2 a; + a;^ -|- 2 a;^ + 2 a;^ 
-I- 4 a; • (5 + a; - 3 a;2 + a;^ 5 a;^) • i^i • 

+ 2x^-{2 -lOx - Ylx^ -lOa;^ + 2x^) ■ FI ■ Fq], (75) 



where Fq and Fi are respectively 2Fi([l/2, 1/2], [1]; a;^) and 2Fi([l/2, 3/2], [2]; x^) 
(closely related, up to a 7r/2 factor, to the usual complete elliptic integrals). Again the 
occurrence of (very simple) elliptic integrals is underlined. Note that L^l^ contributes 
to the solutions of £^y^ in the block L^l^ ■ V^"^ ■ which has the local exponents 

a; = 0, P = 0^ 1^, 2, 

a; = 00, P = 1^, 2^, 3, 

a; = 1, p = -3,-2, 0^ (76) 

x = -\, p=o^2^ 

1 -1- X -I- a;2 = 0, p = 0,l,2,3,4, 5/2. 

There are two solutions analytic at a; = with exponents 1 and 2. 

J)-. 2. On the order-six linear differential operator L^^^ : "Special Geometry" 

Let us write the formal solutions Fg^^ at a; = 0, where the notation [a;^] means that 
the series begins as x^ {const. + •••). There is one set of five solutions and one extra 
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solution analytical at a; = (i.e. four solutions with a log, and two solutions analytical 
at a; = 0): 

5*1 = [x^] ln(x)4 + [x^] \n{xf + [x^] \n{xf + [x"] \n{x) + [a;"], 



^2 - [a;^] ln(:E)3 + [x^] \n{x)^ + [x^] \n{x) + [x% 

S3 = [x^] ln(x)2 + [x^] ln(x) + [x^], (77) 

5*4 = [x^] ln(a;) + [x'^], S5 = [x^], and: 

^6 = 

In view of this structure, the linear differential operator L^^^ cannot be (homomorphic 
to) a symmetric fifth power of the linear differential operator corresponding to the 
elliptic integral. 

(5) 

The next step is to see whether the exterior square of Lq has a rational solution, 

(5) 

which means that ig corresponds to "Special Geometry". With the six solutions 
(|77p . seen as series obtained mod. primes, one can easily built the general solution of 
Ext'iLf) as 

^^'=,^■('5^5 -^p Mp=1,---,6, (78) 

which should not contain log's, fixing then some of the coefficients d^.p. 
For a rational solution of Ext^{Lf^) to exist, the form (free of log's) 

D{x)-Y^d,,,-{Su'^ ~S,^l (79) 

should be a polynomial, where the denominator D{x) reads 

D{x) = • (a; + 1)"^ • (a; - 1)"^' • (1 + a; + x^T^ • (1 + x + a;^ + a;^ + a;^)"^ 

the order of magnitude of the exponents rij being obtained from the local exponents 
of the singularities. With series of length 700, we have found no rational solution for 
Ext\Lf\ 

(51 

Even if ig is an irreducible operator of even order, we have looked for a rational 
solution for its symmetric square. The general solution of Sym?'{L^Q^) is built from 
((771) as 

Y.fk,p-SkSp, fc>p=l, •••,6, (80) 

k,p 

and the same calculations are performed. With some 300 terms, we actually found 
that Sym?{L'^^^) has a rational solution of the fornjj (with Pig6(a:) a polynomial of 
degree 196): 

-PiQeix) 

(a; + l)i0-(x-l)i4.(l + a; + a;2)2i .(l + x + a;2 + a;3+a;4)9' ^ ' 

thus showing that ig does correspond to " Special Geometry" . 

Note that the occurrence (|77|1 of two analytic solutions at a; = 0, for L'f\ which 
is irreducible, is a situation we have encountered in Ising integrals |14[ [T7] . The order 
twelve differential operator (called l!'^2* in [17]) has four analytical solutions at a; = 
and it has been demonstrated that it is irreducible |17j . 

f Note that this form occurs, for a non minimal representative of Lg , in the factorization l|7ip . On 
this point, see the details around (43), (44) in |14| . 
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5. Singular behavior of x'd\x) 

Now that we have obtained all the solutions of the linear differential equations of 
xl^\x) and x'"f\t), analytic at the origin, we turn to the exact computation of their 
singular behavior at the finite singular points. 

To obtain the singular behavior of Xd i^) amounts to calculating the singular 
behavior of each term in (|18p . The details are given in [Appendix D| 

5.1. The behavior of Xd\x) as x ^ 1 

The evaluation of the singular behavior as x — > 1 corresponds to straightforward 
calculations that are given by (IT^ (see (jD.7[) and (ID.Sp ): 

Sol{L^i^){S^ngular,x^l) = -- ^"^ + "i' + (af ^ ~ ( ^°^ + ) 

TT (1 ~ X) V IT 



^ da, . / rs^ fAOn. + a, 

)[bingular, X = ij = 
) / Stt 

'2 ^9r2(5/6)r2(2/3) r2(i/6)r2(i/3)' 



(1 

(3) 

(3) , 57r Stt \ 1 '^2 1 /I \ /Qo\ 

+ " ' ^9r2(5/6)r2(2/3) " r2(i/6)r2(i/3)^J ' +^-Mi--^)- (82) 



When specialized to the combination defining (a^), the singular behavior reads 

'1 Stt 47r \ 1 



X^; {x){Singular,x = 1) = 



3 18r2(5/6)r2(2/3) r2(l/6)r2(l/3)/ l-x 



+ -^.ln(l-x). (83) 

4 TT 

This result agrees with the result determined numerically in Appendix B of |19j . 

One remarks, for the particular combination (|18p giving xf\x), that the most 
divergent term disappears. Note that this is what has been obtained for the 
susceptibility x^'^^ where the singularity (1 — Aw)^'^/'^ of the ODE is not present in 

5.2. The behavior of Xd (^) cls x —1 

The calculations of the singular behavior as a; — >■ —1 rely mostly on connection 
formulae of 2F1 hypergeometric functions, and the results are given below in (jD.9[) 
and (jP.lSp . For the combination (|18p . the singular behavior reads 

x!-!\x){Singular,x=-l) = ^ • ln(l + a;)^ + (^^ _ . ln(l + x), 

which agrees with the result of Appendix B of [TO] , 

5.3. The behavior of xf' {x) as x ^ e±2'^V3 

The result for the singular behavior x'p (x) &s x xq = e='=2^V3 reads: 
Xi-p [Singular, x = xo) = -^-^ e^/^^ . _ ^„)7/2 

= -0.0957529 ••• e"/i2 • (x - xo)^/^ (84) 

This result is in agreement with the numerical result of Appendix B of fl9! namely 
- V2/3 e'^^/12 . a;o)^/^ with b = 0.203122784 • • • 
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6. Singular behavior of xf\x) 

To obtain the singular behavior of Xd'^ [x) amounts to obtaining the singular behavior 
of each term in (IMl) . 

6.1. Behavior of X'^^t) as t 1 

The calculations of the singular behavior of X'i^\{t) as t — > 1 are displayed in 
[Appendix E[ and read 

X^d\i^){Singular,t = I) = 



7r2(l-t)2 7r2(l-t) 
27r2 1-t 27r2 l-t ^ ' 



To compute the singular behavior of X^^3(i) as t — >■ 1 we need the expression 
of the hypergeometric function 4F3([l/2, 1/2, 1/2, 1/2], [1, 1, l];z) as z 1. This 
hypergeometric function is an example of solution of a Calabi-Yau ODE, and explicit 
computations of its monodromy matrices have been given |48j . 

The differential equation for 4J3([l/2, 1/2, 1/2, 1/2], [1, 1, 1]; z) is Saalschiitzian 
and well-poised (but not very- well-poised) . At z = 1 it has one logarithmic solution 
and three analytic solutions of the form 

oo 

5] c„ • (1 - z)". (86) 

n=0 

The Cn satisfy the fourth order recursion relation 

16 n- (n- 1)2 (n - 2) • c„ - 24 {n - l){n ~ 2)(2n'^ - 6n + 5) • c„_i 

-I- 16 (n - 2)2 (3n2 - 12n -I- 13) • c„_2 - (2n - 5)-* • c„_3 ^ 0, (87) 
where c„ = for n < —1. The vanishing of the coefhcient c„ at n = 0, 1, 2, of c„_i 
at n = 1, 2 and c„_2 aX n ~ 2 guarantees that cq, ci, C2 may be chosen arbitrarily. 

The behavior at z = 1 of 4F3([l/2, 1/2, 1/2, 1/2], [1, 1, 1]; z), which is the solution 
of the ODE that is analytic at z = 0, is given in Theorem 3 of Biihring [49] with the 
parameter 

3 4 

' = -T. = 1' (88) 

(i.e. the Saalschiitzian condition [421 1431 H^ ). For completeness we quote this theorem 
which is valid for all p+iFp{[ai, • • • , ap+i], • • • , bp]; z) when the parameter s of 



IS any intege^ s > 
r(ai) • • • r(ap+i) 



m)---T{bp) 



p+iFp{[0'ij • • • J flp+i]; [bi, ■ ■ ■ ,bp]; z) 



s-l 



n=0 



^/<.(l-z)" +^/>.(l-z)" 

n—s 

oo 

(1 - zY ■ + qn ■ ln(l - z)] • (1 - z)", (89) 



n=0 



I Again we emphasise the role of fc-balanced hypergeometric functions. 
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for |1 — 2;| < 1, — TT < arg(l — z) < n and p — 2, 3, • • • where for < n < s — 1 

r< ^ . -,y« r(ai+n)r(a2 + n)(5-n-l)! js - n)k (p) 

^ ' ' r(ai + s)r{a2 + s)nl ' {a, + s)k {a2 + s)k ' " ' ^ ^ 

for s < n 



j> ^ / -^Nn {ai + s)n-s ia2 + S)n-s (fc-n + g)! ^^(p) 

«! \ {ai+s)k{a2 + s)k' ''^ ' ^ ^ 

and it;„ and q„ are such that 



Wn + qn ■ ln(l -z) = (-1)" • ^ ' \, , 

[s + n)\ nl 

X (y 1 ^ ■ ■ k(l + n - fc) + + s + n) 



k=0 



- ip{ai + s + n) - V^(a2 + s + n) - ln(l - z) ^ , (92) 

where (a)„ = a (a + 1) • • • (a + n — 1) is the Pochhammer's symboL The A^^^ are 
computed recursively in [49j as p — 1 fold sums. In particular 

A2) _ (&2 - asjk {bl - Q3)fc .„„x 



and 



(3) _ (&3 + ^2 - a4 - aa + fc2)fc-fc2 (^i ~ «3)fc-fc2 (^3 - a4)fc2 (^2 - a4)fc2 



4^^ = E 



(fc- fc2)!fc2! 

fc2=0 ^ ^ 

_ (&1 + &3 - as - a4)fc (62 + 63 - as - a4)fe , , 

- fc! ^^^^ 

X 3-P'2([&3 - a3, &3 - a4, -fc], [61 + 63 - aa - 04, 62 + &3 - a3 - 84], 1). 
For use in ((63|) we need to specialize to = 1/2, hj — 1, where 

,(3) ^ (l + fc2)fc-fc2 (l/2);.-fc2 (1/2)L 

^ (fc-fc2)!fc2! 

= fc!- 3^^2([l/2,l/2,-fc], [1,1]; 1), (95) 
and for respectively n = and n > 1 

° - 4 A. (3/2)^ ' ^» - ( 1) („)! 1. (3/2)^ ■ (96) 

We note, in particular, the terms 

A^Q^ = 1, = 3/4, = 41/32. (97) 

Using these specializations in (j89p we compute the terms in ^^-^(i) which diverge as 
t — 1. The term (1 — ■ ln(l — t^) cancels and we are left with 

Xg(t)(5^n,.^a.,t ^ 1) ^ -i . + ^ + ^^-^ . (3/> - 4/2>; 



3 77^ 16 



, -In^. (98) 
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Thus, using ([55]) . and we find that the terms in so/(ig*^), which diverge 
as t — 1, are 

(4)^^4)^ 1 



2 



8 7r2 

_(4) 8(34^) -74^^)) , 16 4^) i3I> -4I> 



37r2 37r2 ' ^ V 1 - f 



Ibai'^'^ ~ 16ai*^ . . 16 . 30^,"^ . 9 , 16 



'--H^,) -^-In^C-^), (99) 

TT^ 1 — r 27r-^ 1 — r 

where the constant 3/j^ — 4/1" reads (with 200 digits): 
3/> - AI> = 

-2.212812128930821923547976814986050021481359293357467766171 
630847360232164854964985815375185842526324049358792616932061 
331297671076950376704358248264961101007730925578212714241825 
5205323181711923135264 • • • (100) 
When specializing to the particular combination (|64p . the singular behavior of 
the integral xli\t) reads 

xf\t){Szngular,t= 1) = ^^^-^ • (^1 - ^[64 + 16 • (3/> - 4/2>)]^ 



7 , 16 1 , 2 16 , , 

+ 16^ •^-137 -le^-'^i^r ^'''^ 

This agreeffi with the result determined numerically in Appendix B of |19| . 

We find again and similarly to Xd (^) tha-t the most divergent term disappears 
for the particular combination giving Xd \t). And here again, this is what has been 
observed [12] for the susceptibility x^^^ at the singularity x = 16w^ — 1 which occurs 
in the ODE as z"'^/^ a^d cancels in the integral x'"*''- 

Remark: It is worth recalling that similar calculations for x*-^-*, also based on 
the evaluation of a connection matrix (see section 9 of [H]), require the evaluation of 
a constant that is actually expressed in terms of C(3): 

1 7 

16 7r3 V 9 6^2 



= TirZ3 •(—-«- o-C(3)), (102) 



when the bulk x'"^^ requires a Clausen constant pTJ that can be written as: 

Chin/S) ^ ^•(3-V^(l, 1/3) +3-V(l, 1/6) -8^2)^ (^03) 

It is quite natural to see if the constant 3/j^ — 4/^ given with 200 digits in 
(|100[) . can also be obtained exactly in terms of known transcendental constants 
(C(3), ■■■), or evaluations of hypergeometric functions that naturally occur in 
connection matrices [TT] (see (jF.2p in [Appendix F[ ). This question is sketched in 
[Appendix F[ 

§ Note that there is an overall factor of 2 between this result and the results given in Appendix B 
of |19| which comes from a multiplicative factor of 2 in the series (around t = 0) of xJj^'C*) used 
in | 19| . This applies also to the result of the singular behavior at t = — 1 . 
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6.2. Behavior of Xd*'*(0 ^ ^ ~1 

When t —1 the only singular terms come from Xd^3(0- Furthermore 

the operator of ([5T|) is non-singular at f = —1. Therefore, the only 

singularities in X'f' [t) come from the terms with ln(l — t^) in the expansion 

of 4i^3([l/2, 1/2, 1/2, 1/2], [1,1,1]; at t -1. Thus, from §^ we find that the 

singular part of x'd\t) at t = — 1 reads 

1 -. oo 

xi:L,W= -lx^Ln,(t) - -iln(l-t^). A3 (1-t^r+i), (104) 



n=Q 



with Qn obtained from (|92|) as 



(3/2)^ >^ (-n)fe .3) 



+ (3/2), 

where A'^j^^ is given by (|M| . We know from the exponents of l'^'' at t — — 1 that the 
result has the form (t-l-l)^-ln(t-l-l). Therefore to obtain this term in a straightforward 
way we need to expand the coefficient of ln(l — t"^) to order (1 -|- t)^ in order that the 
term from (l + t) ■ _Df be of order {1 +ty . This is tedious by hand but is easily done 
on Maple and we find that the leading singularity in xli '' (t) at t = — 1 is 

x';,^\Szngular,t^-l) ^ -1- . {I + t)' ■ \n{l + t) , (106) 
which agrees with Appendix B of |T5] . 



7. Conclusion: is the Ising model "modularity" reducible to selected 
(q+i)Fq hypergeometric functions ? 

In this paper we have derived the exact analytic expressions for Xd"^ i^) and x'f^ (t) and 
from them have computed the behavior at all singular points. We have also obtained 
some additional exact results for xf \x) (see section d]). This completes the program 
initiated in [12] where the singularities were studied by means of formal solutions 
found on Maple and numerical studies of the connection problem [llj . In this sense 
we have a complete solution to the problem. However, in another sense, there are still 
most interesting open questions. 

In section 16.11 we used the solution of the hypergeometric connection problem 
|49| which gave the connection constants and as multiple sums. However 
there are special cases, as mentioned in |52| . where it is known by indirect means 
that the series can be simplified, but for which a direct simplification of the series 
has not been found. One example is given by the computation in section 15.21 of the 
singularity of xf'' i^) at < = 1 which we accomplished by means of the reduction (IT5)) 
of a 3F2 function to a product of 2^1 functions. This produced the gamma function 
evaluation of the singularity at x ~ 1 of (|D.8[) . This singularity could also have been 
computed directly from the 3F2 function in (fT6| by use of the Biihring formula (|89)) 
but a reduction of the sums for the required /„ to the gamma function form is lacking. 
There are two suggestions that such a reduction may exist for x'f^t) at i = 1. The 
first is that, by analogy with the corresponding calculation for x (t) ttie bulk [3], 
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the amplitude could be evaluated in terms of C(3)- The second is that evaluations 
of Calabi-Yau gS] hypergeometric functions like 4F3([l/2, 1/2, 1/2, 1/2], [1, 1, 1], z) 
take place. The larger question, of course, is how much the structure seen in x)i ^^id 
■^(") for n = 1, 2, 3, 4 can be expected to generalize to higher values of n. It is the 
opinion of the authors that there is a great deal of mathematical structure of deep 
significance remaining to be discovered. 

These new exact results for the diagonal susceptibility of the Ising model confirm 
that the linear differential operators that emerge in the study of these Ising rt-fold 
integrals, are not only "Derived From Geometry" |13| . but actually correspond to 
"Special Geometries" (they are homomorphic to their adjoints, which means ^27j that 
their differential Galois groups are "special" , their symmetric square, exterior square 
has rational function solutions, ...). More specifically, when we are able to get the exact 
expressions of these linear differential operators, we find out that they are associated 
with elliptic function theory (elliptic functions [53] or modular forms), and, in more 
complicated cases, Calabi-Yau ODEs [HI [53]. This totally confirms what we already 
saw [13] on x^^^ and x^^^. We see in particular, with x'd\^)^ emergence of a 
remarkable order-six operator which is such that its symmetric square has a rational 
solution. 

Let us recall that it is, generically, extremely difficult to see that a linear 
differential operator corresponding to a Calabi-Yau ODE [HI [23] , is homomorphic to a 
q+iFq hypergeometric linear differential operator up to an algebraic pullback. Worse, 
it is not impossible that many of the Calabi-Yau ODEs are actually reducible (up to 
operator equivalence) to q+iFq hypergeometric functions up to algebraic pullbacks 
that have not been found yet. Let us assume that this is not the case, and that 
the Calabi-Yau world is not reducible to the hypergeometric world (up to involved 
algebraic pullback), we still have to see if the "Special Geometry" operators that 
occur for the Ising model, are "hypergeometric" ones, reducing, in fact systematically 
to (selected fc-balanced) q+iFq hypergeometric functions, or correspond to the more 
general solutions of Calabi-Yau equations. 



Acknovifledgment This work was supported in part by the National Science 
Foundation grant PIIY-0969739. We thank D. Bertrand, A. Bostan, J. Morgan, A. 
Okunkov, and J-A. Weil for fruitful discussions. M. van Hoeij is supported by NSF 
grant 1017880. This work has been performed without any support of the ANR, the 
ERG or the MAE. 

Appendix A. Miscellaneous comments on the modular curve (|27p 

Let us introduce other rational expressions, similar to ([T7|) and ([26)) : 

21x^-{l + x) 27x-{\ + xY 
^^("^ = (x + 2)6 ' = ix-lY ' 

where recalling the expression of (|26p one has (for instance): 

OiW = Oi(j) = 0.(--4^). <33(i) = Q'{-YT^) 
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Remarkably the elimination of x between the Hauptmodul Q = Q{x) and 
Q2 = Q2{x) (or Q ~ Q{x) and Qj, = Qz{x)) also gives the same modular curve 

We also have remarkable identity on the same hypergeometric function with these 
new Hauptmodul pullbacks (jA.ip : 

{x+2)- 2F,([\,\]Al];Q{x)) 

1 1, 



2 . (1 + X + a;2)i/2 . ([-, -], [1]; Q^ix)) , (A.l) 



and: 



(l-.x).2^^i([-, -], [1];Q(. 

= (l + :r + :r2)i/2.2Fi([i i], [l];Q3(x)). (A.2) 
The well-known fundamental modular curve |16] 

5%^^^ - 12 • S^w^w^ -{u + v) +375 uv (16 + 16w^ - 4027 um) 
-64:{v + u) ■ {v'^ + U87VU +u^) +2^'^-3^-uv = 0, (A.3) 

corresponding to the elimination of the variable x between the previous Hauptmodul 
(|17p and another HauptmoduI§l QLix): 

(1 +x)'^-x 
{x^ - 14 a; + 1 

should not be confused with the (modular) curve [16] ([27| . 

The new modular curve p7p also has a rational parametrization, (w, v) = 
{Ql(x), Qi(x)), between this last new Hauptmodul (|A.4p and a new simple 
Hauptmodul: 

( 1 +xf -x^ 

[l -X 



i(x) = 108- ^\7'^V ■ (A.5) 



Appendix B. Solution of A^4 analytical at x = 

The solution of A^4 (see ([SO]) ). analytical at x = 0, reads: 

S'o;(x) = (1 - a;)^/^ • p{x) ■ S{x), (B.l) 
where S{x) reads: 

Zi ■ 4F3([i i i i], [1, 1, 1]; P{x)) +Z2 ■ 4F3([| I I [2, 2, 2]; P(a:) 



^3 • 4F3( [| I I [3, 3, 3]; P{x)) + Z4 • ^F^i^-, ^, ^, [4, 4, 4]; P(a;) 



with 



Zi = -512 -^i, Z2 = 128-5^, Z3 = -54-5^, Z4 = -625-5^, 

Ctl 02 «3 «4 

Related by a Landen transformation on a;^/^ see [16) . 
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and 

m = (7x^ - 56x2 + ^^2x - 64) • (1 - xf/"^ 

+ (x - 1) [x^ - 24a;2 + 80 a; - 64), 
712 = (2352 x^ - 472 x^ - 3904 x + 2048 + 19 x^) ■ (1 - xY'"^ 

+ x^ - 125 a;* + 1288 x^ - 4048 + 4928 x - 2048, 
ria = (x^ - 28080 x^ ~ 355 - 52992 x + 17920 + 5750 x^ + 57760 x^) • (1 - xf^"^ 

- 2 (x - 1) (20 x^ - 855 x^ + 6736 x^ - 18992 x^ + 22016 x - 8960), 
714 = (x^ -8x +8)(x'' -64x^ +320x2 - 512x + 256)- (1 -x)^/2 

- 4 (x - 1) (x - 2) (3 X - 4) (x - 4) (x^ - 16 x + 16), 



and 



and 



di = (1 -x) -x^ • ((x- 2) • (x^ - 16x + 16) • (1 -x)i/2 

-2(x- l)(3x-4) (x-4)), 
d2 = (1 -x) -x* • (4- (x- 2) • (1 -x)i/2 + x^ - 8x +8), 
da = (1-x) -x^ • (2(x- 1) - (x-2) • (1 -x)i/2)^ 

d4 — (1 — x) • X®, 



/ , \ 1/2 

p{x) = ((2-x).(l -x)i/2 +2.(x -1)) , 
and where P{x) denotes the puUback ([55)) : 

X^ X 

This solution has the integrahty property [3 2) . Changing x into 64 x the series 
expansion of the previous sohition (jB.ll) has integer coefficients: 

5*0/(64 x) = 128 + 2560X + 116736x2 + 6072320 x^ + 335104000 x'* 
+ 19117744128 x^ + 1114027622400 x^ + 65874638708736 x^ 
+ 3937277209282560 X* + ••• 

(5) 

Appendix C. The linear differential operator C\-{ in exact arithmetic 

The factors occurring in the differential operator read 



^^^-^-^lK (i+V-(i^I!--.- )' ^^-'^ 



^ 2 dx V(l + x)2 • (1 -x)6 -x/' ^ ^ 

Lf = d4 +P1.DI +Pl.Dl+P^-D^ +^, (C.4) 

P4 P4 PA PA 
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with: 

Pi = x^- (1 + a; + a;2) • (a; + 1)^ • (a; - 1)'' • (l60 + 3148 x + 24988 a;^ + 86008 x^ 

+141698 a;^ + 69707 a;^ - 141750 a;^ - 358707 a;^ - 356606 a;* - 1071 x^ + 347302 a;^^ 
+510214 a;^^ + 347302 a;^^ - 1071 a;^^ - 356606 a;^^ - 358707 a;^^ - 141750 a;^'' 
+69707 a;^^ + 141698 a;^* + 86008 x^^ + 24988 x'^° + 3148 ar^^ + 160 a;^^^ , 

P3 = 2x2- (x + l)2- (^-880 - 16620 a;- 126586 .T^- 421558 a;^- 520547 a;'' 

+733378 + 3794648 x^ + 6252130 + 3922367 x^ - 4349032 x^ - 12817741 x^° 
-12881692 x" - 2612141 + 10986996 x^^ + 16830947 x^'' + 12283572 x^^ 
+729267 x^^ - 8919176 x^^ - 10905121 x^* - 5398478 x^^ + 866024 x^° 
+3665682 x^^ + 3069821 x^^ + 1351818 x^^ + 323590 x^* + 36308 x^^ + 1680 x^'^) , 

P2 = 2x-{x-lf ■ (2400 + 38692 x + 228422 x^ + 366806 x^ - 1591741 x^ - 8948446 x^ 

-18137183 - 10301088 x^ + 31576074 x^ + 82978356 x^ + 80098415 x^^ 
-8308172 x" - 123518048 x^^ - 158759046 x^^ - 65285821 x" + 78248130 x^^ 
+152708392x1^ + 124727752x1^ + 26488355x1** - 65301174x1^ - 90679899x2° 
-47527872 x^i + 4032496 x^^ + 27473954 x^^ + 23107094 x^* + 9927812 x^^ 
+2288564 x'^^ + 245416 x^^ + 10800 x^*) , 

Pi = 2 (x-1)- (^-1440 - 15176 X- 3552x2 + 632252 x^ + 3988986 x^ + 11012538 x^ 

+10122851x^-31358640x^- 125311964 X*- 166380144 x^ + 20063039x^° 

+375202188 x" + 523233277 xi^ + 189830162 xi^ - 422078559 x" - 747281488 xi^ 
-440223099 xi^ + 161161298 xi^ + 530901457 xi^ + 491902752 xi^ + 168466049 x'^° 
-168274188 x^i - 282329480 x^^ - 158906808 x^^ - 754525 x^* + 72189798 x^^ 
+61435092 x^s + 25677392 x^^ + 5672988 x^* + 577984 x^^ + 24000 x^°) , 

Po = -3600 - 52880 X- 324108x2- 1147996x3- 1575180 x^ + 8228874 x^ 



= -3600 - 52880 X- 324108x2- 1147996x3- 1575180x^ + 8228874 
+52977905 x^ + 108476130 x^- 739178 x^ - 371064711 x'' - 563202298 a. 
-29824206 x" + 842725375 xi2 + 1075242362 xi^ + 273493047 x" - 909934423 xi^ 
-1189246308 xi'^ - 414338515x1^ + 420114304x1* + 702981552x1^ + 447865799x2 
+30467322x21 - 270639170x22 - 233990685 x23 - 67035676x2^ _^ 45089100 x25 
+61580064 x2<5 + 29851532 x2^ + 7030080 x2* + 714400 x'^^ + 28800 x^^. 



Appendix D. Analysis of the singular behavior of X^'^'(x) 

Let us give some detailed analysis of the singular behavior of xflix) and xfli^) 
around the three singularities: x = +1,-1, 6^2^*73 _ 
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Appendix D.l. Limit of the connection matrix. 

The hypergeometric operator I?^ + ((1 + a + 6) z — c) + ah has two solutions, 
Ml — 2Fi{[a, 6], [c]; z) and U2 = z^^'' • 2-F'i([a + 1 — c, 6 + 1 — c], [2 — c]; z) and 
they connect to z = 1 (see, for instance, (1) on page 108 of [51]) according to the 
connection matrix vahd for c 7^ 1 

2Fi([a,6],[a + 6-c+l]; 1 -2) 
(^_^)c-a-fc. 2F^([c-a,c-6], [c-a-6+1]; 1 -z) 



Ml 




M2 





(D.l) 



with 



C 



Cii C12 
C21 C22 



Cii 
C21 



T{c)r{c - a - b) 
r(c- a)r(c- 6)' 
r(2-c)r(c-a-6) 



(^12 — 

C22 = 



where: 

T{c)T{a + h- c) 

r(a)r(6) ' 

r(2-c)r(a + 6-c) 



(D.2) 
(D.3) 



r(i-a)r(i-6) ' r(a + 1 - c)r(6 + 1 - c)' 

If we now take the limit c — > 1 , which is the case of interest in our problem, the 
connection matrix becomes singular. In this limit, solutions with ln(z) occur, and we 
write M2 = Ml + (1 — c) -112, where mi — 2Fi{[a,h], [1]; z) and 



M2 



Inz • 2-Fi([a, b], [1]; z) 



d_ 

dc 



2Fi([a + l-c,5+l-c], [2-c]; z) 



yielding for the connection of the solutions mi and M2: 

2Fi{[a,b],[a + b]; 1 - z) 



Ml 




U2 





(1-z) 



1-a-f) 



,Fi([l-a,l-6],[2-a-6]; 1 - z) 



(D.4) 



(D.5) 



with 



C21 
C22 



= lim 



C = 

C2I ^ C*!! 



Cii C12 
C2I C22 



where: 



c^l 1 — c 



r(l - a - 6) 
r(l-a)r(l-6) 



(2^(1) -^(1- a) -^(1-6)), 



lim 



c^l 1 - C 

where V(^) = r'(z)/r(z). 



C22 - C12 ^ Tja + h-l) 

r(a)r(5) 



(2VX1) - ^(a) - ^(6)), 



(D.6) 



Appendix D.2. The behavior as x ^ 1 

To evaluate xli-li^) fo^' x — ;> 1 we use the evaluation of 2-F'i([l/2, ±1/2], [1]; x^) for 
a; — >■ 1, and find, from p^ . that its singular part reads 

2 11 

Xrf^^(5mgM?ar,a; = 1) = — + — ln(l - a;). (D.7) 

7r(l — .t)^ 7r(l — x) 2n 

To evaluate xi^al^^) as a; — >■ 1 we specialise (|D.5|) to a = 1/6, 6 = 1/3, 
z = Q, where Q is defined by ([T71) . Then as a; — ^ 1 one has (1 — Q)"^/^ — > 
2/\/3/(l — x) — l/\/3, and, thus, one deduces the singular part of X^'^3(a;) using 
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(fTO|l. or ratheii x^j^^^ix) = Mi(2i^i([l/3, 1/6], [1]; Q)^), where the hnear differential 
operator Mi is defined by (1^ : 



X%iSingular,x = 1) = ^ • — ^ 



(1-.T) 



1 

TT 



(1 -a;)2 



Stt 



9r2(5/6)r2(2/3) r2(i/6)r2(i/3) 



(D.; 



Appendix D.3. The behavior as x —1 

When x — )■ —1 it is straightforward from to obtain 



xl^^i^ingular, x - 



^.ln(l + x) 



(D.9) 



To evaluate xj^s note, when a; — 1, that Q vanishes as Q ~ ^ (1 + x)2. 
However, we cannot directly set Q = in or because we must analytically 
connect the solution analytic at a; = to the proper solution at a; = —1. To do 
this we use the results of |Appendix D.l Using the fact that there is no singularity at 

x = —1/2, we see that we must choose near x — —1/2 



(1-0) 



1/2 _ 



(1 -x) • (1 +2x) • (2 +x) 
2- (l + x + a;2)3/2 



(D.IO) 



which is positive for — 1/2 < a; < and negative for — 1 < a: < —1/2. Therefore, 
for —1 < a; < —1/2, we see that 

r(i/2) 



Ul 



2Fi([l/6,l/3], [1]; Q) 



r(5/6)r(2/3) 



,i^i([l/6,l/3],[l/2];l -Q) 



r(i/6)r(i/3) 

Furthermore for — 1 < x < — 1/2 



2Fi([l/6,l/3], [1], Q) 
V3 



V6)ram ' ' 2-P^i([5/6,2/3], [3/2]; 1 - Q). (D.ll) 

r(-l/2) 



2^ ' 2 ^r(5/6)r(2/3) 
d 



■ C22 + 



r(i/6)r(i/3) 



C2I • Ul 



2^ 2Fi([l/6,l/3], [1]; Q) InQ - — 2Fi([7/6-c,4/3-c], [2-c]; Q)|,^i 



1 



r(i/2) 



■C: 



r(-i/2) 

2 Vr(5/6)r(2/3) r(i/6)r(i/3) 

By using relations (ID.6P we note that 

r(i/2) ~ , r(-i/2) 



r(5/6)r(2/3) 

and thus as a; — )• — 1 



• C22 + 



2Fi([l/6,l/3], [1]; g) ^ ^ (ln(l + a:) -21n2), 

TT 



C21 ] ■ 2^^i([l/6,l/3], [1]; Q). (D.12) 

(D.13) 
(D.14) 



^/3 



,,,,,, C21 = -— • (61n3 + 41n2), 
r(l/6)r(l/3) 27r ^ 



V3 



Again using we find, as x — > —1, that 

3 



X% (Singular, X 



-1 



27r2 



ln2(l + a;) +3 ^^"^„ ^ ■ln(l + x). (D.15) 



t Using l|24ll and the fact that 2i^l([l/3, 1/6], [1]; Q) is a solution of X2. 
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Appendix D.4- The behavior as x ^ g±2iri/3 

When X — g±27"/3 tJ^gn Q — > oo and xi^l becomes singular. Thus to extract this 
singularity we have to connect that solution analytic at a; = to the singularity at a; = 
g±27ri/3^ To do this it is convenient to notice that Q is symmetric about x = —1/2. 
This is seen by letting a; = -1/2 + iy, to obtain Q(y) = {1 + Ay^)^ / {1 - ^y^)^ , 
and defining z by 

z= (l-Q(y))^ = (3^^ _ ^2)3/2 ■ 

Furthermore, as y goes from to V3/2, Q{y) goes from 1 to oo. In the previous 
section we have already connected the solution analytic at x = with the solution 
analytic at a:: = —1/2. 

We rewrite the solutions using (jD.16l) as 

2Fi([l/6, 1/3], [1]; Q) = r{5/6)r(2/3) ' ^Z^^' t^/^^' 

' ^^"^/^^ z- 2Fi([5/6,2/3], [3/2]; z2), (d.17) 



r(i/6)r(i/3) 



^^([7/6,4/3],[2];Q) - 18- ^ 2Fi([l/6, 1/3], [I]; Q). 

These solutions must be connected from y — to y ~ along the straight line 

path X = ^1/2 + iy. On this path is on the negative real axis and, hence, we 
may use the connection formula (2) on page 109 of [51] 

2F,{[a,b],[c];z') - w. iw^ " 1 ■ i-^T" ■ 2fi([«, 1 - c + a], [1 - 5 + a]; z"^) 
1 [0)1 (c — a) 

+ r!1?r~M - • 2F^i[b,l-c + bl[l-a + b]■,z-^). (D.18) 

1 (a)l (c — o) 

Setting z = i z with z real and nonnegative we obtain 

2Fi([l/6,l/3], [1]; Q) = ^(^/3-z)i^0 .£-1/3 . ^^^([l/e, 2/3], [5/6]; -z-^) 

+ 1) • • 2Fi([l/3,5/6], [7/6]; -z''), (D.19) 

and a similar relation for 2-F'i([7/6, 4/3], [2]; Q). Now we note that Q = 1 — z^ = 
1 + z^, and thus (see (fTg])) 

2Fi([l/6,l/3], [1]; +^2i^i([l/6,l/3], [1]; Q) • 2^1 ([7/6, 4/3], [2]; Q) 

can be rewritten in terms of 2^1 hypergeometric functions of argument — this 
last expression going as — 54a/3 • i z'^^/n/SS when z ^ 00. Noting, as a; e^'^^^^, 
that 

we find that the leading singularity at a;o = ^^^^/^ jj-^ -^(3) jg 

x'iliSingular, x ^ Xo) = • e^/^^ . _ ^^^7/2^ ^ 31) 
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Appendix E. Analysis of the singular behavior of Xd-2(^) sls t 1 

To get the singular behavior of X^^lit) as t — >■ 1, we use (12) of page 110 of [5T] 

2Fi([l/2,-l/2],[l];t) ^ (In-IL + ... (E.l) 

and 

,r , , 1 r 1 X 1/16 l-t 16 \ , , 

2Fi([l/2,l/2],[l];t) _. (^ln-_ -2j +••• (E.2) 

Using these in ([T3| we find the result quoted in the text in (|85l) . 

Appendix F. Towards an exact expression for 3/j^ — 4/^ 

The constant 3/{* — 4/^ occurs for x'd\'t) through Xd^l{t) which is ob- 
tained in (j63l) by the action of the linear differential operator ^3 on 
4F3([l/2, 1/2, 1/2, 1/2], [1, 1, 1]; t^). 

The constant 3/{* — 4/^ can then be deduced from the 4x4 connection matrix 
for 4F3([l/2, 1/2, 1/2, 1/2], [1, 1, 1]; t'^). The line of the connection matrix relating 
the solutions at < = to the solutions at t = 1 is 

[A44, -1/2- A44 +2/7r2, A4,3-2i/7r, ^4,4+^/^], 

and the constant 3/j^ — 4/;^ reads 

^•(3/,>-4/>) ^ -i + ^-Ai (F.l) 

The entry A4.1 of the connection matrix is actually the evaluation of the 
hypergeometric function ()59p at t = 1: 

^4,1 = -2.4F3([i, i, i, i], [1, 1, l];l). (F.2) 

There is, at first sight, a "ln(2)" coming from the terms in the Biihring formula [H] 
involving the ip function. This is the "same" ln(2) which appears in the connection 
formulas for E{k) and K(k). However, numerically, these ln(2) contributions in A4 3 
and A4_4 read respectively (a = -1.9453040783 • • • , 7 = 0.5274495683 • • • ): 

^4^3= a + 2 • /3 • ln(2), ^4,4- 7-/3-ln(2), (F.3) 

fi = 0.101321183 ••• 

The fact that these two entries occur through the linear combination ^4^3 + 2 • y44_4 
actually cancels a ln(2) contribution in the expression of the constant ZI^ — 4/^. 

Similar constants (see (|102l) for the bulk (|103p for the bulk x*-'^^) can be 

deduced from entries of the connection matrices (occurring in the exact calculation of 
the differential Galois group [H]), such entries being often closely related to evaluation, 
at selected singular points, of the holonomic solutions we are looking at. When 
hypergeometric functions like (|59p pop out, it is not a surprise to have entries that can 
be simply expressed as these hypergeometric functions at a; = 1 (see (jF.2[) ). Along 
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this line, it is worth recahing that ^(3) (or C(5), ■•■) can be simply expressed in terms 
of a simple evaluation at x = 1 of a q+iFq hypergeometric function |5D] (see also jM]): 

C(3) - 4^3([1,1,1,1],[2,2,2]; 1), (F.4) 

C(5) - 3Y • 6^5([2' 2' 2' 2' 2' ^' ^2' 2' 2' 2' 2^' 

It is thus quite natural to ask if the sums in /{* and can be evaluated in 
terms of known constants such as ^(3) or evaluations (for instance at t = 1) of 
hypergeometric functions. 
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